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ABSTRACT 

Electronic  Oscillators  are  studied  as  an  application  of  the  general 
method  in  the  parameter  plane  of  analyzing  non-linear  systems  with  self- 
excited  oscillations.   Computational  methods  are  developed  which  avoid 
the  lengthy  calculations  required  by  the  general  method  for  the  second 
and  third  order  systems  which  arise.  The  single  parameter  and  multi- 
parameter sensitivity  of  the  frequency  characteristics  are  determined 
and  it  is  demonstrated  that  these  sensitivity  measures  may  be  used  to  de- 
termine the  best  configurations  of  networks  to  be  used  to  minimize  sensi- 
tivity to  specific  parameters.  The  method  employed  suggests  application 
to  the  frequency  characteristics  of  passive  and  active  filters. 

A  specific  problem  illustrating  the  failure  of  the  £  =  0  contour  in 
the  parameter  plane  to  completely  delineate  regions  of  stable  and  un- 
stable solutions  to  the  system  differential  equation  is  included  in  the 
appendix.   No  solution  is  given  but  suggestions  as  to  how  this  apparent 
paradox  arises,  are  made. 


TABLE   OF  CONTENTS 


Monterey,  ca  93943-5101 


Section 
1. 
2. 
3. 


Introduction 

The  Parameter  Plane  Method 

Electronic  Oscillators 

Application  of  the  Parameter  Plane  Method 
Oscillators  as  Feedback  Systems 

4.  The  Sensitivity  Problem  in  Oscillator  Design 

Sensitivity  Definition 

Single  Parameter  Sensitivity 
Multiparameter  Sensitivity 

Frequency  Sensitivity 

Second  Order  System 
Third  Order  System 

Amplitude  Sensitivity 

Sensitivity  in  Specifications 

5.  Analysis  of  Common  RC  Oscillators 

Phase  Shift  Type 

Low  Pass  Configuration 

High  Pass  Configuration 
Other  RC  Oscillator  Types 

Wien  Network 

Bridged  MT"  Network  -  Type  1 

Bridged  "T"  Network  -  Type  2 

Twin  "T"  Network 
Stability  Analysis  of  Limit  Cycles 
Example 


Page 
9 
11 
16 
16 
20 
23 
23 
24 
24 
26 
28 
28 
29 
30 
32 
33 
33 
35 
36 
36 
38 
40 
40 
42 
43 


TABLE  OF  CONTENTS  (cont.) 

Section  Page 

General  Observations  44 

Multiparameter  Sensitivity  45 

Algebraic  Determination  of  Parameters  48 

6.  Analysis  of  Common  LC  Oscillators  51 

Tuned  Plate  Oscillator  51 

Tuned  Grid  Oscillator  53 

Colpitts  Oscillator  53 

Hartley  Oscillator  53 

Sensitivity  Analysis  55 

Example  57 

7.  Conclusions  58 

8.  Bibliography  60 

9.  Appendix  A   A  Basic  Problem  Arising  From  the  Sta-  63 

bility  Analysis  of  Limit  Cycles  In 
the  Parameter  Plane 

Introduction  63 

The  Particular  Problem  63 

Conclusion  65 

Suggestions  69 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Coordinates  of  Point  in  s-Plane  13 

2.  Analysis  of  Stability  of  Limit  Cycles  18 

3.  Block  Diagram  of  Oscillator  21 

4.  Components  of  Typical  £  =  0  Contour  27 

5.  Phase  Shift  Oscillator  -  Low  Pass  Configuration  34 

6.  Phase  Shift  Oscillator  -  High  Pass  Configuration  34 

7.  Wien  Network  37 

8.  Bridged  "T"  Network  -  Type  1  37 

9.  Bridged  "T"  Network  -  Type  2  39 

10.  Twin  "T"  Network  39 

11.  RC  Oscillators  -  Parameter  Plane  41 

12.  RC  Oscillators  -  Frequency  Sensitivity  41 

13.  Tuned  Plate  Oscillator  50 

14.  Tuned  Plate  Oscillator  -  Equivalent  Circuit  50 

15.  Tuned  Grid  Oscillator  52 

16.  Tuned  Grid  Oscillator  -  Equivalent  Circuit  52 

17.  Colpitts  Oscillator  54 

18.  Colpitts  Oscillator  -  Equivalent  Circuit  54 

19.  Hartley  Oscillator  56 

20.  Hartley  Oscillator  -  Equivalent  Circuit  56 

21.  C  *  °  Contour  for  Bridged  "T"  Network  Oscillator  68 


LIST  OF  TABLES 

Table  Page 

1.  Determination  of  Multiparameter  Sensitivity  47 

2.  Parameter  Plane  Analysis  -  Example  66 

3.  Routh-Stability  Criterion  -  Example  67 


ACKNOWLEDGEMENTS 

The  author  is  indebted  to  Dr.  G.  J.  Thaler  of  the  Naval  Postgraduate 
School  for  his  guidance  and  direction  provided  during  the  preparation  of 
this  thesis. 

Sincere  thanks  are  also  extended  to  Dr.  S.  R.  Parker,  also  of  the 
Naval  Postgraduate  School  for  his  helpful  discussions  and  for  his  sug- 
gested references. 


SENSITIVITY  ANALYSIS  OF  SELF-EXCITED 

OSCILLATIONS  IN  NON- LINEAR  SYSTEMS. 

AN  APPLICATION  TO  ELECTRONIC  OSCILLATORS. 

1.    Introduction. 

The  behavior  of  electronic  oscillators  has  been  described  both 
qualitatively  and  quantitatively  in  many  ways.   However,  as  new  models 
and  methods  in  analysis  are  developed  for  the  solution  of  new  problems, 
it  is  worthwhile  to  test  their  application  to  older  problems  to  see  if 
a  simplification  or  unification  of  theory  exists.   It  is  in  this  spirit 
that  electronic  oscillators  are  studied  herein.  They  are  considered  as 
a  class  of  non-linear  active  networks  in  which  self-excited  oscillations 
arise. 

The  principle  analytical  tool  to  be  used  is  the  Parameter  Plane 
Method.  This  application  to  electronic  oscillators  is,  therefore,  a 
practical  application  of  general  theory  developed  by  Siljak[20]  and  Sil- 
jak  and  Stoic[21]  for  self-excited  oscillations  in  non-linear  systems. 

In  practical  oscillator  design  and  application,  the  fundamental 
problem  is  amplitude  and  frequency  stability  with  respect  to  changes  in 
system  parameters.   Second  to  this  may  be  the  accurate  control  of  fre- 
quency by  deliberate  variation  of  system  parameters.  This  problem  is 
analyzed  by  one  parameter  and  multiparameter  sensitivity  analysis.   In 
this  respect,  the  sensitivity  problem  is  approached  in  a  historical  se- 
quence.  It  is  first  considered  as  an  afterthought  to  oscillator  design 
and  then  possible  design  techniques  are  considered  for  use  when  the  sensi- 
tivity is  given  as  an  a  priori  design  specification. 

Several  well  known  oscillator  circuits  of  both  RC  and  LC  type  are 
analyzed  and  sample  problems  are  worked  to  show  how  the  parameter  plane 


representation  and  sensitivity  analysis  can  be  used  to  choose  one  top« 
ology  over  another  to  satisfy  design  specifications. 
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2.   The  Parameter  Plane  Method. 

Before  the  specific  Application  of  the  parameter  plane  method  is 
made,  it  is  worthwhile  to  trace  the  development  of  the  method  as  an 
algebraic  technique  in  control  system  design  and  to  examine  the  basic 
theory. 

Mitrovic[18]  is  generally  credited  with  initiating  the  method  in 
the  early  1950' s.   His  work  was  directed  to  finding  the  roots  of  a  poly- 
nomial as  functions  of  the  coefficients  of  the  first  and  zeroth  order 
terms.   (The  role  of  finding  the  roots  of  a  polynomial  is  a  control 
engineering  oriented  problem,  to  say  the  least!)  He  found  it  possible 
to  map  contours  of  constant  damping  factor,  £,  and  undamped  natural  fre- 
quency, ua  ,  from  the  complex  frequency  domain  into  a  plane  with  the  co- 
efficients as  ordinate  and  abcissa.   However,  considerable  extension  of 
theory  was  required  in  order  to  develop  the  parametric  equations  of 
suitable  contours  in  terms  of  realistic  system  parameters  which  enter 
into  many  coefficients  of  a  characteristic  polynomial  of  a  system.   The 
method,  as  a  graphical  tool  for  finding  roots  and  determining  system 
stability  has  much  in  common  with  the  root  locus  technique,  but  with  the 
facility  to  observe  the  movement  of  roots  simultaneously  varying  with 
respect  to  two  parameters  rather  than  one.   Attempts  have  been  made  to 
extend  the  concept  into  a  parameter  space  to  consider  variation  of  roots 
with  respect  to  more  than  two  parameters,  but  interpretation  of  results 
rests  on  developing  a  greater  facility  with  the  parameter  plane  and  the 
algebraic  expressions  associated  with  it. 

In  the  complex  frequency  domain  (s-plane),  a  root  position  can  be 
characterized  by  a  set  of  cartesian  or  polar  coordinates.  Theory  of 
the  complex  variable  suggests  the  real  and  imaginary  part  of  the  com- 
plex variable,  s,  as  suitable  cartesian  coordinates  and  the  practical 
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criteria  of  damping  factor  and  undamped  natural  frequency  borrowed  from 
second  order  system  analysis  are  suitable,  if  unconventional,  polar  co- 
ordinates.  Figure  1  shows  the  relation  of  these  coordinates. 
It  is  seen  that 

s  =  -  or  +  juj 

or        s  =  -  ud  cos8  +  ju)  sin8 
n        n 

but       £  =  cos6 


hence     sin6  m  /l-g8 
from  which  we  have 

s  a  -  a   +  ju)  (1) 

s  .  -  u)nC  +  jiuyi-j^  (2) 

The  general  nth  order  characteristic  equation  of  a  system  can  be 
represented  by 

n     k 

E  Vs  °  (3) 

k.0  * 

Mitrovic  directed  his  effort  to  developing  expressions  for  the  paramet- 
ric equations  of  contours  of  constant  g  and  u)  in  the  coefficient  plane, 
that  is,  the  chosen  parameters  were  the  coefficients  of  the  first  and 
zeroth  order  terms  in  (3).  The  application  to  finding  the  roots  of 
polynomials  was  direct,  and  in  particular  led  to  a  single  normalized 
graphical  relation  for  all  of  the  roots  of  any  third  degree  polynomial. 
However,  in  application  to  control  systems,  it  was  seen  that  using  co- 
efficients as  parameters  was  less  than  satisfactory  as  a  general  method 
since  the  real  system  parameters  appeared  in  combination  in  many  co- 
efficients. 

Siljak[19]  generalized  the  method  of  Mitrovic  to  the  case  where  the 
real  parameters  appeared  in  linear  combination  in  any  coefficient.   If 
a  and  0  are  two  parameters,  in  general,  we  might  have 
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Coordinates  of  Point  in  s-Plane 
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Subsequent  generalization  by  Hollister[23]  and  others  has  led  to  solu- 
tions for  systems  in  which  the  parameters  appear  as  products  or  to  the 
second  power.  To  date,  this  has  been  sufficient  to  allow  the  applica- 
tion of  the  method  to  a  large  number  of  problems  in  control  system  and 
network  analysis. 

At  the  heart  of  the  method  is  the  simple  fact  that,  since  s  is  a 
complex  variable,  then  any  equation  f(s)  s  0  has  a  real  and  imaginary 
part  which  are  independently  zero.  This  gives  rise  to  two  simultaneous 
equations  in  any  pair  of  parameters  a   and  $  in  terms  of  the  coordinates 
of  the  s-plane,  <C  and  %  from  (2)  or  a  and  u>  from  (1)). 

Recent  work  by  Siljak[20]  has  shown  that  this  essentially  linear 
technique  can  be  extended  to  yield  accurate  results  in  the  prediction  of 
limit  cycles  and  the  stability  of  limit  cycles  in  non-linear  systems. 
The  frequency  of  a  limit  cycle  can  be  determined  and,  in  some  cases,  the 
transient  response  of  a  non-linear  system  may  be  constructed  from  the 
parameter  plane  representation.  The  application  is  subject  to  limita- 
tions at  the  present  time  but  is  based  on  the  conjecture  that,  if  a  lin- 
ear system  can  be  represented  by  the  collection  of  stationary  roots  of 
the  characteristic  polynomial,  then  a  non-linear  system  may  be  character- 
ized by  a  set  of  moving  roots.   Knowledge  of  the  dynamics  of  root  behav- 
ior should  then  yield  knowledge  of  the  system  behavior  in  the  time  do- 
main.  Of  course,  no  definition  exists  of  a  "dynamic  root"  at  this  point 
in  time,  since  the  whole  concept  of  a  characteristic  polynomial  of  a 
system  is  based  on  the  essentially  linear  analysis  technique  of  the  La- 
place transform.   Yet,  the  notion  of  moving  roots  is  producing  meaning- 
ful results  under  certain  restrictions. 
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The  application  of  the  parameter  plane  method  to  non-linear  systems 
is  just  such  an  extension.  Whereas,  in  a  linear  system  we  speak  of  a 
particular  root  system  being  characterized  by  an  "M-point"  which  is  fix- 
ed in  the  parameter  plane  by  the  choice  of  a  parameter  pair,  we  must  now 
speak  of  a  "dynamic  M-point",  the  locus  of  the  values  of  the  parameters 
as  dictated  by  the  nature  of  the  non-linearity.   This  is  the  technique 
to  be  used  in  the  analysis  of  electronic  oscillators. 

From  Mitrovic's  original  idea  has  been  developed  a  comprehensive 
analysis  tool.   From  time  to  time  one  hears  objections  to  the  method  as 
being  unnecessarily  complicated  or  even  redundant  in  the  face  of  other 
analytic  methods.  Granted,  the  method  applied  in  general  is  best  used 
where  a  general  purpose  computer  is  available  to  plot  the  parameter 
plane  contours,  yet,  for  many  problems  there  are  simplifications  which 
show  the  utility  of  the  method  while  avoiding  the  complexity.  The  appli- 
cation of  the  method  to  analysis  of  the  sensitivity  problem  in  frequency 
of  oscillators  is  just  such  a  case. 
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3.    Electronic  Oscillators. 

In  oscillators  producing  a  nearly  sinusoidal  wave-form,  amplitude 
stabilization  is  achieved  by  the  onset  of  some  non-linear  property  of 
an  element  of  the  network.  This  could  be  the  effect  of  saturation  or 
the  change  in  gain  of  the  active  element,  and  in  some  practical  oscilla- 
tors it  has  been  brought  about  by  using  a  lamp  as  a  non-linear  resist- 
ance to  affect  the  bias  of  a  variable  transconductance  type  vacuum  tube 
to  control  the  gain.   In  the  usual  analysis,  it  is  assumed  that  the 
oscillation  can  be  described  in  terms  of  the  location  of  a  pair  of  poles 
of  the  closed  loop  transfer  function  on  the  imaginary  axis  of  the  com- 
plex frequency  plane.   Such  a  delicate  balance  is  highly  unlikely  to 
last  very  long  in  a  truly  linear  system.   The  application  of  the  para- 
meter plane  analysis  removes  the  necessity  of  the  assumption,  provides 
a  new  insight  into  the  stability  or  instability  of  the  oscillation 
which  is  properly  a  limit  cycle,  and  then  provides  a  validation  for  the 
above  assumption  for  the  case  of  oscillations  with  a  nearly  sinusoidal 
wave- form. 

Application  of  the  Parameter  Plane  Method 
In  this  analysis,  it  is  assumed  that  there  is  only  one  non-linear 
element.   One  of  the  two  variable  parameters  is  assigned  to  this  non- 
linear property.   If  there  are  two  non-linearities,  a  variable  para- 
meter designation  is  assigned  to  each.   If  there  are  more  than  two  such 
non-linearities,  then  assumptions  must  be  made  as  to  a  representative 
single-value  to  be  used  for  all  but  two  non-linearities  or  extensions 
must  be  made  into  a  "parameter  space"  where  one  would  consider  a  "dy- 
namic M-surface"  rather  than  the  dynamic  M- locus  considered  in  the  two- 
dimensional  parameter  plane.   In  the  case  of  electronic  oscillators, 
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there  is  usually  only  one  significant  non-linearity  in  the  system,  and 
this  will  be  the  case  in  all  of  the  examples  studied  in  this  thesis. 
The  second  variable  parameter  is  assigned  to  a  system  parameter  of  in- 
terest and  all  other  parameters  are  fixed  at  some  reasonable  value.  By 
mapping  the  stable  and  unstable  regions  of  the  parameter  plane,  and  by 
superimposing  the  excursions  of  the  parameter  assigned  to  the  non-lin- 
earity, it  is  possible  to  determine  the  existence  of  a  limit  cycle,  to 
predict  its  stability  and  to  predict  the  frequency. 

In  the  s-plane,  the  stable  and  unstable  regions  are  bounded  by  the 
8  s  j«)  or  £  at  0  line.   This  contour  can  very  easily  be  mapped  into  the 
parameter  plane  without  recourse  to  the  extensive  general  methods.   By 
substituting  s  =  juu  into  the  characteristic  equation  of  the  oscillating 
system  ,  the  equation  can  be  readily  divided  into  a  real  and  imaginary 
part,  which  equated  separately  to  zero  yield  two  simultaneous  equations 
in  the  parameters,  say  a  and  g ,  as  functions  of  the  frequency  oo.   These 
are  the  parametric  equations  of  the  contour  £  m   0  or  a  =  0  in  the  para- 
meter plane. 

Since  there  is  only  one  non-linearity  in  consideration,  the  M-locus 
is  a  straight  line.   The  intersection  of  this  line  with  the  £  a  0  con- 
tour indicates  the  existence  of  a  limit  cycle  at  the  frequency  given  by 
the  parameter  uu  at  the  point  of  intersection.   A  qualitative  argument 

It  is  at  this  juncture  that  it  is  quite  proper  to  object  from  a 
mathematical  standpoint  that  this  whole  structure  violates  some  very 
basic  concepts.   A  characteristic  equation  arises  from  the  Laplace  Trans- 
formation of  the  system  differential  equation.   If  any  coefficient  or 
coefficients  are  time-varying  (as  they  most  assuredly  will  be  in  the  case 
considered),  they  cannot  be  removed  from  the  transformation  integral  and 
it  is  not  correct  to  carry  such  coefficients  into  a  "characteristic 
equation".   I  cannot  answer  such  an  objection,  indeed  anyone  who  can  has 
a  good  handhold  on  the  problem  of  analyzing  non-linear  systems  as  systems 
with  moving  roots.   Yet  this  is  the  basis  on  which  I  will  proceed  and  let 
the  results  themselves  speak  for  the  validity  of  the  assumption  and  beg 
for  a  more  complete  mathematical  justification! 


17 


M-locus 


Stable 


Unstable 


Stable 


c-  ° 


a 


Figure  2 
Analysis  of  Stability  of  Limit  Cycles 
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may  then  be  given,  based  upon  the  nature  of  the  non-linearity  as  to 
whether  the  limit  cycle  is  stable  or  not. 

Figure  2  shows  a  possible  £  =  0  contour  for  a  system  in  which  the 
gain  of  a  non-linear  active  device  is  a  function  of  the  amplitude  of 
the  signal.   In  this  case,  as  the  amplitude  increases  the  gain  decreases. 
The  plane  is  divided  into  three  regions  for  positive  real  values  of  the 
parameters  a   and  @ .   The  point  X  on  the  end  of  the  M-locus  is  determined 
by  the  maximum  possible  gain  and  the  point  Y  by  the  minimum  possible  gain. 
If  the  oscillation  should  start  with  initial  conditions  placing  the  M- 
point  in  the  region  C,  the  amplitude  will  decrease,  since  the  oscillation 
is  stable,  thus  causing  the  gain  to  increase  and  the  M- point  moves  toward 
the  £  =  0  contour  at  a.   If  the  initial  conditions  are  such  as  to  put  the 
M-point  in  the  region  B,  the  oscillation  is  unstable  and  will  grow  in 
amplitude,  causing  the  gain  to  decrease  and  the  M-point  to  move  down  the 
locus  to  the  intersection  with  the  £  =  0  curve  at  a.   Thus,  regardless 
of  whether  the  oscillation  starts  with  M  in  region  A  or  B,  the  amplitude 
will  change  until  a  steady  state  condition  of  a  stable  limit  cycle  is 
produced  with  amplitude  and  frequency  which  characterize  the  point  of  in- 
tersection at  a.   Should  the  amplitude  be  such  that  the  M-point  is  in  re- 
gion C,  the  region  is  stable,  the  amplitude  will  decrease,  the  gain  will 
increase  and  the  M-point  will  move  up  toward  X  and  away  from  the  inter- 
section at  b.  The  intersection  at  b  is  also  a  limit  cycle,  but  it  is 
clear  that  no  oscillation  can  ever  be  established  at  that  frequency  in 
the  system  described  since  the  tendency  at  b  is  either  to  have  the  os- 
cillation collapse  or  move  to  the  stable  limit  cycle  conditions  of  b. 

It  is  seen  from  the  above  analysis,  that  the  steady  state  conditions 
are  such  that  the  self-excited  oscillation  may  be  characterized  by  the 
intersection  of  the  M-locus  with  the  £  ■  0  contour.   But  this,  in  itself 
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is  equivalent  to  assuming  that  an  oscillator  may  be  analyzed  by  consid- 
ering its  oscillation  to  be  characterized  by  a  pair  of  poles  on  the  juu 
axis  of  the  s-plane.   Hence  the  application  of  the  parameter  plane  meth- 
od justifies  the  assumption  made  in  earlier  analyses  and  provides,  in 
terms  of  non-linear  theory  an  insight  into  the  operation  of  the  oscil- 
lator.  It  should  be  noted  that  the  earlier  analyses  were  restricted  to 
the  steady  state  oscillation  whereas  the  parameter  plane  theory  has  suf- 
ficient information  to  describe  the  transient  period  after  turning-on 
during  which  the  oscillation  is  moving  to  the  stable  limit  cycle.   If 
the  duration  of  this  transient  or  the  excursions  from  steady  state 
brought  about  by  changes  in  system  parameters  are  of  interest,  expres- 
sions for  these  changes  have  been  developed. [20] 

Oscillators  as  Feedback  Systems 
In  general,  any  oscillator  that  provides  continuous  self-regulation 
can  be  considered  as  a  feedback  system.   The  diagram  in  Figure  3  shows 
such  a  general  system  in  which  the  forward  element  is  considered  to  be 
an  active  element  in  which  the  gain  is  a  function  of  the  amplitude  of 
the  signal.   The  general  characteristic  equation  of  such  a  system  from 
block  diagram  algebra  is 

1  -  0H(s)  *  0  (5) 

From  this  simple  relationship  the  Barkhausen  criteria  for  oscillation 
are  immediately  apparent,  i.e. 

|gH(s)|  =  1  (6) 

/BH(s)  =  0  (7) 

The  loop  gain  must  be  unity  and  the  total  phase  shift  around  the  loop 
must  be  some  integral  multiple  of  2rr . 
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Figure  3 
.lock.  Diagram  of  Oscillator 


21 


In  many  practical  oscillators,  the  active  element  is  a  single  stage 
vacuum  tube  or  transistor  which  provides  an  inherent  n  phase  shift.   In 
such  a  case  the  feedback  network  will  provide  an  odd  integral  multiple 
of  tt  phase  shift. 

The  feedback  network  itself  may  take  on  many  forms  and  topologies 
and  there  is  a  fundamental  problem  in  oscillator  design  to  choose  a  suit- 
able feedback  network  topology  and  to  choose  parameters  to  suit  specifi- 
cations.  The  parameter  plane  method  of  analysis  and  comparison  of  pro- 
perties of  different  topologies  and  a  sensitivity  analysis  of  different 
topologies  provides  a  sound  basis  for  choice. 

It  is  convenient  to  consider  oscillators  to  be  divided  into  two 
classes,  those  incorporating  resonant  circuits  as  the  frequency  deter- 
mining feedback  networks  and  those  incorporating  only  resistive  and 
capacitive  elements  in  the  feedback  circuit.   The  latter  are  important 
in  the  development  of  integrated  circuits  where  an  inductive  element 
cannot  be  manufactured  in  the  circuit  with  the  ease  that  a  resistive  or 
capacitive  element  can.   In  this  thesis,  RC  type  oscillators  are  con- 
sidered in  general  with  an  idealized  active  element  considered  to  be  a 
saturating  amplifier  with  infinite  input  impedance  and  zero  output  im- 
pedance, or  an  amplifier  with  deadzone.   The  LC  oscillators  considered 
are  analyzed  on  the  basis  of  the  Norton  equivalent  circuit  of  a  vacuum 
tube  active  device.   Extension  to  transistorized  circuits  is  a  matter 
only  of  substituting  a  suitable  model  and  redetermining  the  character- 
istic equations  of  the  oscillators.   The  analysis  is  exactly  the  same 
as  that  to  be  developed. 
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4.    The  Sensitivity  Problem  in  Oscillator  Design. 

In  practical  use,  oscillators  are  required  to  produce  and  maintain 
a  waveform,  often  sinusoidal,  of  fixed  frequency  and  magnitude.   If  the 
application  requires  that  the  frequency  and/or  amplitude  be  variable, 
the  variation  must  be  accurately  controlled.   However,  there  are  many 
parameters  in  an  oscillating  system  which  are  subject  to  variation  in 
a  non-deterministic  manner.   The  reduction  of  sensitivity  of  amplitude 
and  frequency  to  these  uncontrolled  changes  is  a  basic  consideration. 
The  sensitivity  analysis  of  oscillators  from  the  parametric  equations 
from  the  parameter  plane  method  provides  a  ready  comparison  of  topolo- 
gies and  the  effect  of  parameter  values  on  sensitivity  to  parameter 
changes.   In  the  following  analysis,  the  effect  of  variation  of  para- 
meters is  the  primary  consideration.   It  is  recognized  that  the  stabil- 
ity of  the  oscillator  signal  also  depends  on  phase  shifts  in  the  ampli- 
fiers and  changes  in  harmonic  content  of  the  signal.   In  many  practical 
circuits,  additional  feedback  compensation  is  made  to  make  the  ampli- 
tude relatively  insensitive  to  frequency  changes.   This  adaptation  is 
not  considered.   However,  it  will  be  noted  that  under  some  circumstances 
certain  unadapted  systems  have  an  inherently  low  amplitude  sensitivity 
to  frequency  variation. 

Sensitivity  Definition 
Sensitivity  is  a  measure  of  the  change  in  an  operating  character- 
istic or  representation  of  a  system  with  respect  to  the  variation  of  a 
system  parameter  or  parameters.   The  classical  sensitivity  problem  deals 
with  small  parameter  changes  and  the  sensitivity  operator  is  differen- 
tial in  nature.   Bode!"l]  is  responsible  for  the  first  sensitivity  defin- 
ition.  Latterly,  several  definitions  have  been  given  and  each  has  its 
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own  utility  in  the  context  of  suitable  problems. 

Single  Parameter  Sensitivity 

The  basic  single  parameter  definition  of  sensitivity  is 

T   change  in  T 
x   change  in  x 

The  change  in  either  T  or  x  can  be  a  relative  change,  i.e.  a  normalized 
or  percentage  change,  or  an  absolute  change.   In  the  measure  of  changes 
in  parameter  values,  when  the  parameters  are  the  values  of  lumped  net- 
work elements  with  percentage  tolerances,  to  speak  of  a  relative  change 
is  useful.   In  the  definition  of  a  sensitivity  function  which  is  the 
measure  of  the  movement  of  a  pole  or  a  zero  of  a  transfer  function,  the 
absolute  value  of  the  change  is  of  significance.   In  the  case  of  measur- 
ing the  sensitivity  of  frequency  of  oscillators  to  changes  in  parameters 
of  the  oscillating  system,  the  relative  change  of  both  frequency  and 
parameter  is  of  interest.  As  the  problem  is  developed,  it  may  be  deter- 
mined that  a  particular  sensitivity  definition  results  in  a  simplifica- 
tion of  expression  or  calculation  which  favors  one  definition  or  another 
In  this  thesis,  the  Horowitzf4]  definition  is  used  as  follows: 


x  =  dx7x"  ■  T  *  dx"  =  d  In  x  K*} 


This  definition  will  be  shown  to  lead  to  some  simplifications  in  the 
calculation  of  frequency  sensitivity  of  oscillators.  The  inherent  sin- 
gularity at  the  null  value  of  the  parameter  should  be  noted,  but  it 
causes  no  difficulty  if  it  is  kept  in  mind  that  this  definition  is  a 
measure  of  a  percentage  change  in  T  due  to  a  percentage  change  in  the 
parameter  x. 

Multiparameter  Sensitivity 

In  many  cases,  it  is  highly  likely  that  no  single  parameter  will 
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cause  all  of  the  changes  in  the  function  of  interest.   A  single  measure 
of  the  cumulative  effect  of  the  variation  of  several  parameters  is  re- 
quired.  Goldstein  and  Kuof 15]  extended  the  single  parameter  definition 
of  sensitivity  into  the  multiparameter  case.   The  total  logarithmic  dif- 
ferential can  be  written: 

j  i      m   £  5  In  T   .  ,  ,rtV 

d  InTr  .  ff*-s d  In  x.  (9) 

i  o  In  x.       i 

By  considering  the  set  of  fractional  parameter  increments  as  a  vector, 
dY  m  (d  In  x.,  d  In  x2,  — ,  d  In  x  ) ,  and  setting  y.  =  In  x.,  the  above 
total  differential  can  be  written 

d  In  T  =  V[ln  T(yv  yfi   -.,  yn)]«  dY  (10) 

If  we  fit  this  expression  to  our  single  parameter  definition,  the  multi- 

-T 

parameter  sensitivity,  S  ,  is  given  by 

S*  a  V  In  T  (11) 

The  magnitude  of  this  sensitivity  vector  can  be  expressed  by  the  norm. 

Still  this  multiparameter  definition  leaves  something  to  be  desired 
for  a  general  definition.   In  the  manufacture  of  integrated  circuits, 
for  instance,  all  of  the  components  are  laid  down  at  the  same  time  under 
identical  physical  conditions.  The  value  of  the  components  can  only  be 
realized  approximately,  but,  since  they  are  all  laid  down  together  a 
high  degree  of  correlation  is  expected.   In  this  case,  it  might  well  be 
that  the  sensitivity  with  respect  to  some  parameter  is  of  opposite  sign 
with  respect  to  another.   Indeed,  good  design  would  take  maximum  advan- 
tage of  such  a  situation  to  reduce  the  variance  in  overall  performance 
criteria  for  the  circuit.   In  this  case,  it  is  not  reasonable  to  take  a 
sensitivity  measure  which  adds  the  squares  of  the  individual  sensitiv- 
ities as  is  the  case  where  the  norm  of  the  sensitivity  vector  is  taken 
as  the  multiparameter  sensitivity.   Leef26]  and  Biswas[31]  have  recog- 
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nized  this  point  and  in  Lee's  concept  of  "Sensitivity  Groups"  and  in 
Biswas*  inclusion  of  correlation  terms  in  the  definition,  a  more  rea- 
sonable sensitivity  measure  for  the  multiparameter  case  is  obtained. 
In  the  case  of  networks  composed  of  individual  components  rather  than 
an  integrated  circuit,  the  designer  might  safely  assume  independence  of 
the  parameter  variations  for  which  the  Goldstein  and  Kuo  multiparameter 
sensitivity  definition  is  adequate. 

Frequency  Sensitivity 
The  parameter  plane  representation  of  an  oscillator  has  been  dis- 
cussed.  With  a  definition  of  single  parameter  sensitivity  in  hand,  it 
is  now  possible  to  see  how  the  concept  of  frequency  sensitivity  may  be 
expressed  graphically  in  the  parameter  plane  and  algebraically  from  the 
parametric  equations  and  from  the  characteristic  equation.   Figure  4 
shows  a  typical  £  s  0  contour  from  a  parameter  plane  drawn  for  an  oscil- 
lator.  From  the  previous  analysis,  it  was  seen  that  the  intersection 
of  the  M-locus  of  the  non-linearity  gives  the  frequency  of  the  limit 
cycle.   The  amplitude  of  the  limit  cycle  may  be  given  if  a  known  rela- 
tionship exists  between  the  amplitude  and  the  defined  non-linear  para- 
meter.  It  is  seen  that  every  point  on  the  £  =  0  contour  is  a  potential 
limit  cycle  point  and  the  frequency  of  the  limit  cycle  is  given  implic- 
itly by  the  parametric  equation  of  the  contour.   Since  the  contour  is 
described  by  the  parametric  equations,  a(u))  and  j3(tu),  it  is  possible  to 
give  the  slope  of  the  graph  as  follows: 

i-jfrte  <») 

from  which  it   is  seen  that 

a  d£      (£  dj|  or  duo  M_* 
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Figure  4 
Components  of  Typical  £  =  0  Contour 
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or 

s^a^s*  (14) 

Hence  it  is  seen  that,  if  B  is  the  parameter  assigned  to  the  variable 
gain  in  the  oscillator,  due  to  the  non-linearity,  the  frequency  sensi- 
tivity with  respect  to  the  variable  parameter,  a,  is  given  algebraically 
from  the  logarithmic  derivative  of  the  parametric  equation  a   =  Qf((w)  or 
from  the  slope  of  the  £  =  0  contour  itself  it  there  is  a  concomitant 
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of  S0. 
m 

Since  the  parametric  equations  of  the  £  m    0  contour  are  ultimately 
determined  from  the  characteristic  equation,  it  will  be  of  some  utility 
to  determine  a  relationship  between  frequency  sensitivity  with  respect 
to  variation  of  the  parameters  and  the  coefficients  of  the  characteristic 
equation. 

Second  Order  System 

For  the  general  second  order  system, 

2 
a_s  +  a.s  +  aQ  ■  0  (15) 

Substituting  s  a  juo,  it  is  seen  that 

U)  a  —  and  a.  a  0   .  (16) 

a2        * 

This  gives  the  characteristic  frequency  and  conditions  for  oscillation 

(two  roots  on  the  ju)  axis).   The  frequency  sensitivity  with  respect  to 

any  parameter  x,  S  ,  is  then  determined  from 

2  In  u)  ■  In  a_  -  In  a?  (17) 

SMa  ^(Sa0  -  Sa2)  (18) 

x   2V  x     x 

Third  Order  System 

For  a  general  third  order  system 
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3      2 
a3s  +  a-s  +  aj^s  +  aQ  a  0  (19) 

Substituting  s  m    juj,  it  is  seen  that 
2   a3   a2 

•  -  -r  =  ;r  (20) 

ai     ao 

from  which  it  follows  that 

tii        1     al  a3  1     a0         a2 

S*  -  i(S        -   S   J)  =  4<S        "  S      )  (21) 

x2xx2Nxx  v      ' 

It  must  be  borne  in  mind  that  the  basis  of  this  definition  in  this  con- 
text is  the  argument  that  there  is  a  self-adjusting  parameter  in  the 
system  causing  the  steady-state  M-point  to  follow  the  £  =  0  curve.   Since 
this  parameter  will  be  included  in  some  or  all  of  the  coefficients  of  the 
characteristic  equation,  the  partial  derivative  of  a  coefficient  with  re- 
spect to  a  parameter  must  also  include  the  partial  derivative  of  the 
adaptive  parameter  in  accordance  with  the  chain  rule. 

Amplitude  Sensitivity 
Since  the  effective  gain  of  the  non- linear  active  element  in  this 
discussion  of  oscillators  is  considered  to  be  a  monotonic  function  of 
amplitude  of  the  signal,  and  since  the  effective  gain  is  specified  by 
the  intersection  of  the  M-locus  with  the  £  ■  0  curve,  it  is  a  temptation 
to  speak  generally  of  the  amplitude  sensitivity.   However,  the  whole 
concept  of  amplitude  with  respect  to  oscillations  in  non-linear  systems 
needs  better  definition.   From  the  chain  rule  it  is  possible  to  show 
that  if  the  effective  gain  is  a  function  of  amplitude  only, 

S3  m   S^.SA  (22) 

A  8 

where  S   is  the  amplitude  sensitivity  with  respect  to  frequency.   S   is 

derived  from  the  parametric  equations  of  the  £  m   0  curve  and  S\  depends 

entirely  on  the  relation  of  effective  gain  to  amplitude.   It  is  in  the 
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relationship  between  amplitude  and  effective  gain  that  there  is  diffi- 
culty.  If  there  is  sufficient  low  pass  filtering  in  the  feedback,  net- 
work, the  describing  function  is  a  reasonable  relation  to  use  for  the 
non-linear  gain.   In  many  cases,  in  practical  oscillators  this  is  not 
the  case.   This  statement  is  equivalent  to  saying  that,  unless  other 
steps  are  taken,  the  waveform  produced  cannot  be  nearly  sinusoidal  but 
will  have  harmonic  distortion  to  some,  possibly  unacceptable  degree.   In 
order  to  reduce  harmonic  distortion  in  such  cases,  it  is  possible  to  ar- 
range that  the  linear  gain  portion  of  the  non-linear  gain  curve  is  just 
high  enough  to  guarantee  that  the  end  of  the  M-locus  lies  just  inside 
the  unstable  region  of  the  parameter  plane.   In  such  cases,  the  ampli- 
tude is  limited  by  the  signal  level  at  which  the  non-linear  property  of 
the  active  device  is  effective  to  reduce  the  loop  gain  to  unity.   This  in 
itself  will  cause  the  amplitude  to  be  frequency  insensitive  but  will  re- 
quire a  degree  of  sophistication  in  the  oscillator  design,  if  the  fre- 
quency is  to  be  variable. 

Sensitivity  in  Specifications 
Sensitivity  specifications  may  be  made  in  many  forms.   The  most 
general  specification  for  a  fixed  frequency  oscillator  is  likely  to  be 
that  the  frequency  shall  not  vary  more  than  a  specified  percentage  due 
to  any  reasonable  cause.   Due  to  the  stochastic  nature  of  the  variation 
of  parameters,  this  necessarily  implies  that  the  tolerance  cannot  be  met 
by  all  oscillators  but  shall  be  met  by  some,  presumably  very  high  number. 
In  addition,  it  is  clear  that  the  specification  is  essentially  a  multi- 
parameter specification  involving  the  calculation  of  sensitivities  due 
to  each  parameter  and  the  probability  density  for  each  parameter.   If, 
however,  there  is  one  predominant  variable  parameter  in  a  system,  until 
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such  time  as  the  effect  of  that  parameter  is  reduced  to  the  order  of  the 
variation  of  other  parameters,  it  is  valid  to  ignore  the  lesser  effects 
and  treat  the  problem  essentially  as  a  single  parameter  sensitivity 
problem. 

After  a  study  has  been  made  of  single  parameter  and  multiparameter 
sensitivity  functions  for  practical  oscillator  circuits,  it  will  be  pos- 
sible to  comment  further  on  suitable  forms  for  stating  sensitivity  speci- 
fications. 
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5.   Analysis  of  Common  RC  Oscillators. 

The  general  problem  considered  in  this  analysis  is  to  find,  on  the 
basis  of  the  comparison  of  parameter  plane  curves  and  sensitivity  func- 
tions, the  best  choice  of  a  feedback  network  consisting  entirely  of  re- 
sistive and  capacitive,  passive  devices  which  will  yield  the  most  stable 
frequency  in  the  sense  of  minimum  frequency  sensitivity.   The  single 
parameter  variation  which  is  likely  to  outweigh  the  others  is  the  load 
resistance.   This  resistance  may  be  the  input  resistance  of  the  active 
device  which  is  likely  to  vary  and  the  input  resistance  of  the  succeed- 
ing stage  in  any  device  in  which  the  oscillator  is  used.   In  actual  fact, 
the  nature  of  the  variable  parameter  considered  is  of  no  importance  to 
the  analysis  since  the  method  is  identical  for  every  parameter.   However, 
conclusions  drawn  in  the  following  examples  are  valid  only  for  the  para- 
meter values  chosen  and  for  the  particular  variable  parameter  considered. 
Each  set  of  parameters  requires  its  own  analysis. 

The  parameter  plane  curves  are  drawn  for  a  normalized  frequency 
X-   RC(D,  thus  making  the  analysis  applicable  to  any  frequency  range. 
Arbitrarily,  to  provide  a  basis  for  comparison  of  different  networks,  all 
resistances  are  chosen  equal  to  R  other  than  the  variable  parameter,  aR, 
and  all  capacitances  equal  to  C.   It  is  recognized  that  this  may  not  lead 
to  the  best  choice  of  parameters  from  other  points  of  view.   Mehta[16] 
has  shown  that  this  choice  is  not  optimum  to  reduce  sensitivity  of  fre- 
quency to  phase  change  in  the  active  element.   However,  this  does  not 
affect  the  method,  of  which  this  material  is  exemplary. 

The  characteristic  equation  of  each  RC  oscillator  was  found  by  the 
method  of  determining  the  transfer  function  of  the  feedback  network  from 
two-port  network  analysis  and  then  substituting  in  1  -  (3H(s)  =  0,  (5), 
the  characteristic  equation  in  general  from  the  feedback  block  diagram. 
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Phase  Shift  Type 
The  phase  shift  oscillator  takes  advantage  of  the  inherent  ft  phase 
shift  of  a  single  or  odd  number  of  stages  of  the  active  element  and 
provides  an  additional  rr  phase  shift  from  the  feedback  network  to  meet 
the  zero  phase  shift  criterion  for  oscillation.  The  feedback  network 
may  be  either  a  high  pass  or  a  low  pass  configuration. 

Low  Pass  Configuration 

Figure  5  shows  the  general  circuit  in  which  the  parameter  of  in- 
terest in  the  comparative  analysis  is  marked  (*) .  The  characteristic 
equation  in  general  is: 

s  (R]R2R3R4C1C2C3)  + 

2 

S  (R.  R_R, C, C~  ">  R-i  R9R/.^i  "o    Ri  R^RoC,  Cn   "r 

R2R3R4C2C3  +  R1R3R4C1C3  +  R1R3R4C2C3)  + 

S(R1R4C1  +  R2R4C2  +  R3R4C3  *  R2R3C2  +  Wl  + 

R1R4C3  +  R2R4C3  +  R1R2C1  +  R1R4C2  +  R1R3C2)  + 

(Rx  +  R2  +  R3)  +  (l-g)R4=  0  (23) 

For  the  particular  example,  setting  R.  =  R_  »  Ro  =  R>  ra  =  aR»  an<* 

C,  as  C?  =  C-  -   C,  and  setting  p  ss  RCs,  we  have  the  reduced  characteristic 

equation, 

ap3  +  (5a  +  l)p2  +  (6a  +  4)p  +  (3  +  a  -  «0)  =  0  (24) 

By  setting  p  ss  jX,  to  plot  the  £  ss  0  contour  in  the  parameter  plane,  we 
have  the  parametric  equations  of  the  contour, 

a  -  WTl  <25) 

_  *  %   nf  ±   * 
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Figure  5 
Phase  Shift  Oscillator  -  Low  Pass  Configuration 
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Figure  6 
Phase  Shift  Oscillator  -  High  Pass  Configuration 
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and  calculating  directly  from  the  characteristic  equation  coefficients 
the  frequency  sensitivity  with  respect  to  a    is  seen  to  be 

S*>  ■ —  (27) 

a  3a  +  2  l  ' 

The  sign  preceding  [3  is  a  measure  of  the  requirement  that  there  be  art 
phase  shift  in  the  active  device  since  the  network  is  introducing  a  n 
phase  shift.   Figure  11  is  a  comparative  plot  of  the  £  *  0  contours  in 
the  phase  plane  for  all  the  examples  of  RC  oscillators  considered.   In 
order  to  compare  the  curves  of  the  phase  shift  types  with  the  others,  the 
carve  is  plotted  for  |{3  |  . 

High  Pass  Configuration 

This  is  the  most  common  phase  shift  oscillator  configuration.  The 
network  is  shown  in  Figure  6.  The  general  characteristic  equation  is: 

B3[(1-P)R1R2R3C1C2C3]  + 

s2(R1R2C1C2  +  R2R3C2C3  4-  R^C^  4- 

R,  R„C-,  C0  T  R-.  R-C,C-j  +  R,  R-C-jC, ,)  T 

s(R1C1  +  R2C2  +  R3C3  +  R1C2  +  R2C3)  +  1  «  0  (28) 

Setting  RjSsRjsR,  R^=aR,  C^C^sC-rrC,  we  obtain  the  reduced  characteristic 
equation  for  the  particular  problem  as  follows 

(a-ap)p3  +  (3f3a)p2  +  (4+a)p  +  1  *  0  (29) 

Once  again  the  third  order  method  is  used  to  obtain  the  frequency  sensi- 
tivity from  the  coefficients  of  the  characteristic  equation.   One  can 
use  the  coefficients  of  the  odd  or  even  powers,  but,  in  this  case,  a 
prior  survey  would  show  that  the  odd  coefficients  are  functions  of  both 
Oi   and  jB  ,  whereas  the  coefficients  of  the  even  terms  are  much  simpler. 
From  the  general  expression, 
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s*  ,  I(sao  -  sa2 

and  the  fact  that  S  "  =  0,  it  is  seen  that 


ASau   -  S°z)  (30) 

a   2V  a     a  v  ' 


a 


S*a  -Isa2=  -■_&. 


(31) 


a  2a  "   2(1  +  a) 

Setting  p  =  j&  the  parametric  equations  of  the  £  a  0  contour  are: 

1  -  "M^ 
o   =   yj  (32) 

-^'r^V*^-1  <33> 

where  once  again  the  sign  before  p  indicates  the  ft  phase  shift  in  the 
active  element  and  the  curve  is  plotted  for)j3|. 

Other  RC  Oscillator  Types 
In  other  types  of  RC  Oscillators,  the  frequency  of  oscillation  is 
determined  by  the  zero  phase  shift  criterion.   Each  of  the  following 
networks  has  a  zero  phase  point  where  the  phase  changes  from  leading  to 
lagging,  or  vice  versa.  The  determination  of  the  characteristic  equation, 
parametric  equations  of  the  £  «  0  contour  and  sensitivity  of  frequency 
to  a  single  parameter  change  proceed  as  before. 

Wien  Network 

The  network  is  shown  in  Figure  7.  The  general  characteristic  equa- 
tion is: 

s2(R1R2C1C2)  +  s[R1C1  +  R2C2  +  (1-0)^^]  +1=0  (34) 

Setting  R,ssR,  Rj»aR,  C^C^aC  we  obtain  the  reduced  characteristic  equa- 
tion 

ap2  +  (1  +  2a  -  ap)p  +1.0  (35) 

from  which,  setting  p  =  jX,  we  find 
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Figure  7 
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Figure  8 
Bridged  "T"  Network  -  Type  1 
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(36) 

(37) 

SJ  =  "  |  (38) 

It  is  worth  noting  that  the  network  is  extremely  simple  but  at  the 
cost  of  having  the  sensitivity  of  frequency  with  respect  to  any  parameter 
change  being  a  constant  value  higher  than  that  in  either  of  the  two  phase 
shift  configurations  thus  far  considered.   For  the  sake  of  comparison, 
the  frequency  sensitivities  are  plotted  in  Figure  12. 

Bridged  "T"  Network  -  Type  1. 

Since  there  is  a  symmetry  between  resistive  and  capacitive  elements 
in  Figures  8  and  9,  the  configuration  of  Figure  8  is  arbitrarily  called 
type  1  for  the  sake  of  identification.   The  general  characteristic  equa- 
tion is: 

s2[(l-B)R1R2R3C1C2]  + 

s[(l-3)(R2R3C1  +  R2R3C2)  +  R1R3C2  +  R]R2C1  +  R!R2C2^  + 
Rx  +  (1-B)R3=  0  (39) 

Substituting  as  before,  the  reduced  characteristic  equation  is 

(a-a8)p2+  (2+3a-2aB)p  +  (l+a-a0)  -  0  (40) 

from  which  it  is  determined  that  the  parametric  equations  of  the  £  =  0 
contour  are 

a  =  ^  (41) 

e  =  *g?  («) 

and  the  frequency  sensitivity  with  respect  to  a  is: 
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Figure   9 
Bridged   "T"  Network   -  Type   2 
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Figure  10 
Twin  "T"  Network 
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a   =  2-k*  (43) 


Bridged  "T"  Network  -  Type  2. 

This  network  is  shown  in  Figure  9.   The  general  characteristic  equa* 
tion  is: 

s2[(l-p)R1R2R3C1C2]  + 

s[(l-p)(R1R3C1  +  R2R3C1)  +  R1R3C2  +  Rj_R2C2]  + 

Rx  +  R2  +  (1-P)R3  =  0  (44) 

Substituting  as  before,  the  reduced  characteristic  equation  is: 

(a-ap)p2  +  (l+3a-2a£)p  +  (2-ky-ap)  =  0  (45) 

from  which  it  is  determined  that  the  parametric  equations  of  the  £  =  0 

contour  are 

3  +'X2 
«  =  ITT?  (46) 

P  =  JTT^  (47) 

*nrf         <£  I*  (48> 

and      a   =  (a+l)(c*-3) 


Twin  "T"  Network 

This  network  is  shown  in  Figure  10.  The  general  characteristic 
equation  is: 

s3[(l-P)R1R2R3R4C1C2C3]  + 

S2r(l-P)(R1R3R4C1C2  +  R2R3R4C1C2)  +  R1R2R4C2C3  +  R^^C^  + 

R,R«R~C„C~  T  R-.R^R,C«C~  ■  Ri  RoR/,^1  ^oJ 

s[(l-g)(R3R4C]L  +  R3R4C2)  +  R2R3C2  +  R2R3C!  +  R1R3C2  +  R1R3C1  + 

R2R4C2  +  R1R4C2  +  R1R4C3  +  R1R2C3^  + 
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[Rx  +  R2(1-P)R4]  =  0  (49) 

Substituting  as  before,  the  reduced  characteristic  equation  is: 

(a-c*e)p3+  (2+5a-2*e)p2  +  (5+5a-2ap)p  +  (2+<y-ap)  =  0  (50) 

from  which  it  is  determined  that  the  parametric  equations  of  the  £  =  0 
contour  are 

a   "   3-3,/  <51> 

2/6  +   9X4-2]X2-10  ,„. 

P  =  2Xb-9/2+   2 (52) 

and  £  •  ±  [S^]"1  (53) 

ca        4*6  +  2*^  +  4X2  ,_,. 

where     S/2  =  -  2/8  +  ^.tfU.tft   +  i  <54> 

It  is  worth  noting  that  this  last  analysis  maybe  about  as  difficult 
an  analysis  as  one  might  consider  making  using  algebraic  operations  with- 
out the  aid  of  a  computer.   If  the  computer  is  available,  the  generalized 
sensitivity  equations  of  Siljak  and  Stoic[2l]  may  be  used. 


Stability  Analysis  of  Limit  Cycles 
As  noted  in  the  theory  of  self-excited  non-linear  oscillations  in 
the  parameter  plane,  the  stable  regions  of  the  parameter  plane  as  parti- 
tioned by  the  £  =  0  contour,  must  be  examined.   If  we  assume  that  the  non- 
linearity  is  such  that  the  parameter,  (3 ,  is  a  monotone  decreasing  func- 
tion of  the  amplitude  of  the  signal,  it  is  quite  easy  to  take  test  points 
above  and  below  the  £  si  0  curve  to  determine  the  stability  of  the  limit 
cycles.   On  this  it  can  be  shown  that  both  phase  shift  networks  and  the 
Wien  Network  give  rise  to  stable  limit  cycles  and  both  bridged  "T"  types 
and  the  twin  "T"  network  give  rise  to  unstable  limit  cycles.   If  a  non- 
linearity  is  chosen  in  which  the  gain  increases  as  the  amplitude  increases, 
such  as  a  dead-zone  type,  the  latter  three  networks  give  rise  to  stable 
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limit  cycles.   Appendix  A  describes  a  basic  problem  in  the  parameter 
plane  method  arising  from  this  analysis 

Example 

Consider  the  details  of  a  possible  design  analysis  in  the  parameter 
plane  based  on  the  curves  developed  thus  far.  An  oscillator  is  to  be 
built  using  an  amplifier  with  a  linear  portion  gain  of  70.  A  phase 
shift,  high  pass  configuration  has  been  chosen.   It  is  desired  to  mini- 
raize  sensitivity  to  the  network  terminal  resistance.   (See  Figure  6.) 
What  should  or  be  and  what  is  the  effect  on  the  sensitivities  with  respect 
to  the  other  parameters  if  ot    is  changed  from  an  arbitrary  reference  at 
a   »  1? 

The  sensitivities  are  calculated  from  the  coefficients  of  the  char- 
acteristic equation.   It  is  seen  that 

S*  =  -  -ka2  (55) 

x     2  x 

where  x  is  any  network  parameter  since 

Sa°  =  0  (56) 

x 

The  sensitivities  at  a  —  1  are: 


S^  =  -  0.41b  S^  s  -  0.333         S^    =   -  0.250 

Ri  R~                   R~ 

rfj  m   -  0.250  s^  =  "  0-333         So  =  "  0.416 

Cl  C2                   C3 


(57) 


The  maximum  possible  sensitivity  with  respect  to  any  single  parameter  is 
0.5  (see  Figure  12).   From  Figure  12,  sensitivity  with  respect  to  R  is 
minimized  by  making  a   as  small  as  possible.   However,  references  to  Fig- 
11  shows  'immediately,  that  3  increases  as  ot    is  decreased,  hence  the 
linear  portion  gain  of  70  is  an  upper  bound  on  (3  ,  implying  a  lower  bound 
on  a  .   A  suitable  value  for  en    is  0.3  for  which  3  a  55.  The  sensitivities 
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at  a   =  0.3  are: 

S*  a  -  0.462 
Rl 

S®  =  -  0.423 
R2 

S*  m   -  0.295 

S^  =  -  0.333 
L2 

S*  =  -  0.115 
R3 


S*  *  -  0.372 
C3 


(58) 


It  should  be  noted  that  although  the  required  minimization  has  been  ef- 
fected within  the  constraints  of  the  problem,  the  sensitivity  with  re- 
spect to  three  others  has  been  increased!  This,  apparently,  is  yet 
another  application  of  the  universal  principle  that  you  can't  get  some- 
thing for  nothing! 

Reference  to  Figures  11  and  12  also  shows  immediately  that,  in  the 
problem  above,  second  thought  should  be  given  to  the  choice  of  a  high 
pass  phase  shift  network.  The  low  pass  phase  shift  network  would  have 
provided  minimization  of  sensitivity  to  R  much  lower  than  the  0.115 
figure  reached  for  any  value  of  a   greater  than  2.5.  The  value  of  R^  is 
very  likely  to  be  determined  by  the  amplifier  specified.   If  it  is  suit- 
able to  choose  R  less  than  R_/a  for  the  required  frequency,  then  a  change 
to  the  other  phase  shift  configuration  is  recommended.   If  this  choice 
is  not  suitable  (the  required  C  may  be  very  large)  then  it  may  well  be 
that  the  phase  shift  high  pass  configuration  may  be  the  only  one  capable 
of  handling  the  low  input  impedance  with  a  low  sensitivity. 

General  Observations 
A  scrutiny  of  Figures  11  and  12  reveals  many  of  the  characteristics 
of  the  RC  Oscillators  studied.   For  those  oscillators  with  a  portion  of 
the  £  a  0  contour  almost  parallel  to  the  a   axis,  it  is  seen  that  dB /d* 
is  very  nearly  zero.   Since  the  frequency  sensitivity  with  respect  to 
a    is  not  zero  in  these  portions,  it  must  be  that  S   is  very  nearly  zero. 
This  implies  that  the  amplitude  sensitivity  is  very  nearly  zero  in  these 


44 


regions,  regardless  of  the  relationship  between  p  and  amplitude.   Hence, 
to  minimize  amplitude  sensitivity  to  change  of  frequency  in  a  variable 
frequency  oscillator,  it  is  worthwhile  to  design  in  a  flat  portion  of 
the  £  =  0  contour. 

It  is  immediately  seen  from  the  parameter  plane  and  sensitivity 
curves  which  configurations  can  meet  an  a  priori  single  parameter  sensi- 
tivity, and  the  required  amplifier  gain  is  also  shown.   The  M-locus  for 
the  non-linear  amplifier  indicates  whether  it  is  feasible  to  consider  a 
configuration,  since  there  may  not  be  an  intersection  with  the  £  =  0 
contour,  or  if  there  is  such  an  intersection,  the  indicated  limit  cycle 
may  not  be  stable,  or  the  sensitivity  specified  may  require  an  impossibly 
high  gain.  A  set  of  curves  can  be  drawn  for  any  parameter  of  interest 
using  suitable  preselected  values  of  the  other  parameters. 

Multiparameter  Sensitivity 
The  comparisons  in  Figure  12  were  based  on  the  single  parameter 
sensitivity  of  frequency  with  respect  to  a  terminal  resistance  in  the 
feedback  network.   Just  so  long  as  the  variation  of  this  terminal  resis- 
tance is  the  predominant  variation  of  parameters  in  the  oscillator,  then 
it  is  valid  to  make  comparisons  of  configurations  based  on  this  single 
sensitivity.   However,  it  is  apparent  that,  when  there  is  a  likelihood 
of  several  parameters  varying  with  an  equal  effect,  then  the  multipara- 
meter sensitivity  measure  must  be  taken.   Provided  the  parameter  varia- 
tions are  uncorrelated,  the  Goldstein  and  Kuo  definition  is  satisfactory. 
Since  this  definition  of  sensitivity  is  a  vector  sensitivity,  the  maxi- 
mum sensitivity  is  given  by  the  norm  of  the  sensitivity  vector.   From  the 
definition 

d  In  T  «  V(ln  T) • (d  In  x. ,  ~,  d  In  x  )  (10) 

1  n 
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rT»  V(ln  T)  «  |   }n  T     x.  +  •••  +  *   }n  T     x  (11) 

where  S  olnx.l  dlnxn  v/ 

1  n 

and  SLx  ■    1^1   ■   tt  <Sx   >^  <59> 

i    i 

On  this  basis,  it  is  possible  to  compare  the  multiparameter  sensitivi- 
ties of  the  two  types  of  phase  shift  networks  considered  in  the  discus- 
sion. 

It  is  assumed  that  the  system  parameters  have  been  chosen  in  ac- 
cordance with  some  other  performance  criteria  and  it  is  desired  to  com- 
pare the  multiparameter  sensitivities  of  the  two  phase  shift  types  and 
Wien  bridge  type  of  feedback  network.  Arbitrarily,  for  the  sake  of  ex- 
ample, it  is  assumed  that  all  resistances  are  equal  to  R  and  all  capac- 
itances are  equal  to  C.  This  does  not  fix  the  frequency,  hence  the 
problem  is  not  completely  abstracted  from  reality!  The  reduced  charac- 
teristic equation  is  determined  for  each  parameter  considered  in  turn  to 
be  variable  or,  to  be  even  more  general,  the  general  single  parameter 
sensitivities  are  determined  from  the  general  characteristic  equation 
for  all  system  parameters.  The  norm  of  the  sensitivity  vector  is  de- 
termined from  the  components  (single  parameter  sensitivities).  Table  1 
shows  the  coefficients  of  the  reduced  characteristic  equations  for  each 
of  the  three  types  considered.   It  should  be  noted  that  this  is  not  a 
minimum  sensitivity  but  a  demonstration  of  the  feasibility  of  getting  a 
multiparameter  sensitivity  index  and  the  fact  that  it  can  be  used  as  a 
basis  of  comparison  of  feedback  networks. 

From  Table  1,  it  can  be  deduced  that  for  a   =  1,  the  multiparameter 
sensitivity  of  the  phase  shift,  low  pass  configuration  network  is  less 
than  either  of  those  of  the  other  two  types  considered.  This  is  similar 
to  the  trend  already  apparent  in  the  single  parameter  sensitivity  to  the 


46 


Table   1 
Determination  of  Multiparameter  Sensitivity 

Phase  Shift   Oscillator-Low  Pass   Configuration 


Parameter 


R, 


R, 


o? 
a 

Of 

a 
a 
a 
a 


2  1  0 

5a  +  1  6a  +  4  a  +  3  -  p 

4a  +  2  4a  +  6  a  +  3  -  B 

4a  +  2  4a+6  a  +  3  -  B 

5a  +  1  6a  +  4  a  +  3  -  ap 

4a  +  2  3a  +   7  4  -  B 

4a  +  2  4a  +   6  4  -  P 

4a  +  2  3a  +   7  4  -  B 


Phase  Shift  Oscillator-High  Pass  Configuration 

Rx  a  -  ap           5a  +   1  2a  +  3 

R_  a  -  aB           4a  +  2  2a  +  3 

R3  a  -  aB            3a  +  3  a  +  4 

C.  a  -  aB            3a  +  3  a  +  4 

C„  a  -  aB           4a  +  2  2a  +  3 

C3  a  -  aB             a  +   1  2a  +^3 

Wien  Network  Oscillator 


For  all   parameters,   S    =   -  -r  .      Hence,   S  =   1.000 
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S  a   0.833 
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terminal  resistance  for  a  a  1  shown  in  Figure  12.   Once  again,  curves 
might  be  drawn  for  comparison  purposes,  but  it  is  recommended  that  a  gen- 
eral  purpose  computer  be  used  if  one  is  available  since  the  calculation 
of  one  point  on  the  curve  represents  considerable  time  spent  in  calcula- 
tions. 

It  should  be  noted  that  several  assumptions  have  been  made  in  de- 
termining this  sensitivity  index  for  many  varying  parameters.   It  has 
been  assumed  that  each  of  the  parameters  is  independent,   a  questionable 
assumption  if  one  is  considering  the  effect  of  ambient  temperature  on 
frequency  stability  since  it  is  highly  likely  that  all  similar  components 
would  suffer  similar  variation  when  subjected  to  the  same  variation  in 
temperature.   It  has  been  assumed  that  the  variation  (relative  or  incre- 
mental) for  each  parameter  is  the  same.   In  the  case  of  the  stochastic 
nature  of  the  variation  of  parameters,  this  assumption  is  equivalent  to 
assuming  that  all  parameters  have  the  same  variance,  valid,  perhaps  for 
the  passive  elements  in  the  network  but  questionable  when  one  considers 
the  terminal  resistance  of  the  network.   If  each  of  the  incremental  var- 
iations of  parameters  is  known,  a  better  assessment  of  the  incremental 
change  in  frequency  to  be  expected  can  be  made.   However,  to  keep  the 
work  of  determining  system  sensitivity  to  a  minimum  it  should  be  noted 
again  that  the  predominance  of  any  single  parameter,  renders  variations 
due  to  other  parameters  secondary  and  the  problem  might  be  considered  a 
single  parameter  sensitivity  problem  with  much  time  saved. 

Algebraic  Determination  of  Parameters 
Thus  far,  the  discussion  has  been  based  on  the  assumption  that  the 
oscillator  design  has  been  substantially  completed  for  several  configur- 
ations and  then  the  sensitivity  analysis  is  applied  to  determine  which 
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design  is  best  suited  to  the  requirements.   In  each  case  it  is  possible 
to  determine  which  configuration  represents  the  minimum  sensitivity  to  a 
parameter  of  interest  but  there  is  no  indication  that  this  sensitivity 
is  the  minimum  for  that  system.   Indeed,  it  might  be  possible  to  find  a 
minimum  sensitivity  set  of  parameters  for  one  of  the  configurations  re- 
jected that  might  make  it  more  suitable.   Lest  it  be  considered  that 
this  is  a  simple  min-max  problem,  note. that  there  are  physical  con- 
straints on  the  solution  of  the  minimum  sensitivity  problem,  such  as 
the  fact  that  all  parameters  must  be  real,  positive  numbers;  capacitors 
larger  than  about  100  mfd  and  resistors  larger  than  22  megohms  are 
specialty  items;  amplifier  gains  for  simple  one  or  two  stage  active  de- 
vices are  seldom  over  100;  the  nature  of  the  sensitivity  functions  al- 
ready analyzed  show  functions  which  do  not  have  minima  for  finite  values 
of  the  argument.   In  some  very  simple  cases,  it  is  possible  to  derive 
general  sensitivity  functions.   By  finding  the  conditions  for  a  minimum, 
one  equation  of  n  possible  equation  is  determined  to  establish  the  values 
of  the  n  parameters.   The  frequency  required  is  not  considered  as  one  of 
the  equations  in  this  context  for  RC  oscillators  since  our  entire  devel- 
opment has  been  based  on  a  normalized  frequency  which  is  determined  by 
the  choice  of  R  and  C  after  the  design  is  optimized. 

One  need  only  look  at  the  nature  of  the  coefficients  of  the  general 
characteristic  equations  in  the  previous  examples  to  realize  that  the 
sensitivity  functions  are  of  such  a  complexity  that  it  would  take  several 
hours  to  set  up  the  min-max  problem.   It  appears  that  an  algebraic  deter- 
mination to  the  choice  of  parameters  for  a  minimum  of  sensitivity  is  a 
problem  well  suited  to  a  computer  solution. 
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Figure   13 
Tuned-Plate   Oscillator 


F 


•g   e 
m  g 


P         <      L 


Figure  14 
Tuned-Plate  Oscillator  -  Equivalent  Circuit 
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6.   Analysis  of  Common  L-C  Oscillators. 

The  elements  of  the  application  of  sensitivity  analysis  to  oscil- 
lator design  have  already  been  developed  and  demonstrated  for  several 
RC  oscillators.   However,  for  the  sake  of  completeness  and  to  point  out 
some  of  the  essential  differences  the  class  of  oscillators  incorporating 
resonant  circuits  is  considered.   The  active  device  in  these  oscillators 
is  the  non-linear  device  once  again  but  in  this  case  there  is  not  the 
same  idealization.   The  active  device  is  considered  to  be  a  vacuum  tube 
and  the  Norton  equivalent  (current  source)  is  used  in  the  analysis.   The 
justification  for  using  the  linear  analysis  on  an  essentially  non-linear 
problem  is  as  before;  the  oscillation  is  a  stable  limit  cycle  and  can  be 
characterized  in  the  parameter  plane  by  the  steady-state  position  of  the 
M-point  moving  on  the  £  =  0  contour.   The  characteristic  equations  aris- 
ing are  second  and  third  order  allowing  the  application  of  the  simplified 
calculation  of  sensitivities  from  the  characteristic  equations. 

Tuned  Plate  Oscillator 

A  resonant  tank  circuit  in  the  plate  circuit  is  inductively  coupled 
to  the  grid  circuit  to  provide  the  necessary  positive  feedback  to  sus- 
tain oscillation.   The  circuit  diagram  and  equivalent  network  are  shown 
in  Figures  13  and  14.   The  general  characteristic  equation  is: 

s2(L1C)  +  s(Liyo  +  Cr  +  gmM)  4-  (1  +  ryQ)  =  0  (60) 

where     yQ  *  -  +  -  (61) 

It  is  worth  noting  that  the  characteristic  equation  coefficients  do  not 
depend  on  L„  but  only  on  the  nqitual  inductance  M.   In  order  for  the  cir- 
cuit  to  oscillate,  the  mutual  inductance,  M,  must  be  negative,  i.e.  there 
is  a  right  and  a  wrong  way  to  hook  up  the  coupling  transformer.   This  is 
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Figure  15 
Tuned-Grid  Oscillator 
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Figure  16 
Tuned-Grid  Oscillator  -  Equivalent  Circuit 
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seen  from  the  coefficient  of  the  first  order  term  which  must  vanish  for 
oscillation  and  cannot  do  so  unless  the  mutual  inductance  is  negative. 
This  might  also  be  stated  that  the  vacuum  tube  introduces  a  ft  phase 
shift  which  must  be  compensated  by  a  ff  shift  through  the  coupling  trans- 
former in  order  to  satisfy  the  unity  loop  gain,  zero  phase  Barkhausen 
criteria  for  oscillation. 

Tuned  Grid  Oscillator 

A  resonant  tank  circuit  in  the  grid  circuit  inductively  coupled  to 
the  coil  in  the  plate  circuit  provides  the  required  feedback.   The  cir- 
cuit diagram  and  equivalent  circuit  are  shown  in  Figures  15  and  16.   The 
general  characteristic  equation  is 

■^CCLjLj  -  M2)yQ]  +  s2[C(Liryo  +  L2>]  + 

s(Cr  +  L,y  +  g  M)  +  1  =  0  (62) 

1  o   Dm 

Colpitts  Oscillator 

This  oscillator  and  the  Hartley  oscillator  are  very  similar  in  that 
the  plate  signal  is  coupled  back  into  the  grid  by  a  balanced  resonant 
circuit.   The  Hartley  is  more  difficult  to  analyze  due  to  the  inevitable 
resistance  of  the  inductors  and  also  due  to  the  mutual  inductance  be- 
tween the  inductors,  which  does  not  exist  in  the  Colpitts  case.   The 
Colpitts  circuit  and  equivalent  circuit  are  shown  in  Figures  17  and  18. 
The  general  characteristic  equation  is: 

s3(LC1C2)  +  s2(LC2yQ  +  rC^)  + 

s(rC2yo  +  Cx   +  C2)  +  (yQ  +  gm)  =  0  (63) 

Hartley  Oscillator 

The  circuit   and  equivalent  diagrams  are  shown  in  Figures  19  and 
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Figure  17 
Colpitts  Oscillator 


J 
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Figure  18 
Colpitts  Oscillator  -  Equivalent  Circuit 
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20.  The  general  characteristic  equation  is: 

s3(CgmLlL2  +  CyoLlL2  -  Cg/  -  Cy/)  + 

s2(CL1r2gm  +  CL2rlgm  +  CL^^  +  Cl^r^  +  C^  +  CL2  +  2CM)  + 

s(Cgmrlr2  +  Cyorlr2  +  Vo  +  Crl  +  Cr2  "  gmM)  +  (r2yo+  l)   =  ° 

(64) 

Sensitivity  Analysis 

At  first  sight,  the  sensitivity  analysis  appears  rather  straight 

forward.   The  equations  are  second  or  third  order  and  the  sensitivities 

with  respect  to  any  parameter  can  be  determined  from  the  coefficients. 

However,  the  use  of  the  Norton  equivalent  of  the  vacuum  tube  is  a  step 

closer  to  reality  in  this  analysis  and  there  are  now  two  parameters 

which  are  not  independent.  To  make  matters  a  little  worse,  these  two 

parameters  are  also  the  non- linear  g  and  r  of  the  vacuum  tube.  A  crit- 

m      p 

ical  analysis  of  the  vacuum  tube  characteristics  will  show  which  of  the 
two  (g  or  r  )  is  the  least  linear  and  presumably  the  most  effective  in 
stabilizing  the  amplitude  of  oscillation.  This  parameter  should  be 
chosen  as  the  non-linearity  in  the  parameter  plane  representation.   In 
the  case  where  the  second  non- linear  parameter  is  a  parameter  of  inter- 
est, there  may  be  a  method  of  constructing  the  M-locus  to  find  the  in- 
tersection with  the  £  *  0  contour,  but  only  if  there  is  a  known  rela- 
tion between  the  parameters  and  the  signal  amplitude.  Little  error  is 
introduced  by  assuming  the  least  effective  non-linear  parameter  to  be 
constant  for  the  sake  of  the  analysis. 

Another  important  difference  presents  itself  in  that  the  RC  net- 
works lent  themselves  to  a  very  convenient  normalization  which  made  a 
single  analysis  applicable  to  oscillators  at  any  frequency.   In  the  case 
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Figure  19 
Hartley  Oscillator 
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Figure  20 
Hartley  Oscillator  -  Equivalent  Circuit 
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of  the  LC  oscillators,  the  characteristic  equations  do  not  lend  them- 
selves to  such  a  simplification.   Even  so,  sensitivity  analysis  can 
lead  to  some  immediate  conclusions  concerning  these  oscillators. 

Example 
The  element  most  likely  to  vary  in  these  oscillators  is  the  induc- 
tance.  If  the  tuned  plate  and  tuned  grid  oscillators  are  compared,  it 
is  seen  that  for  the  tuned  plate  oscillator 

whereas,  for  the  tuned  grid  oscillator, 

If  the  variations  of  L.  and  L_  in  the  tuned  grid  oscillator  are  uncor- 
rected, the  multiparameter  sensitivity  with  respect  to  L,  and  L~  will 
be  less  than  ~  .   However,  if  they  are  correlated  the  multiparameter  sen- 
sitivity will  be  equal  to  —  .  This  is  derived  from  the  Goldstein  and 
Kuo  definition  of  sensitivity.   Hence,  it  appears,  that  there  is  likely 
to  be  a  smaller  frequency  sensitivity  due  to  changes  in  inductance  in 
the  case  of  the  tuned  grid  circuit. 


57 


7.   Conclusion. 

Parameter  plane  analysis  and  sensitivity  analysis  for  the  frequency 
characteristics  of  oscillators  can  be  carried  out  quite  readily  for  sys- 
tems with  second  and  third  order  characteristic  equations,  without  re- 
sort  to  lengthy  and  generalized  parameter  plane  expressions.   Comparison 
of  the  £  =  0  contours  in  the  parameter  plane  gives  a  useful  basis  of  com- 
parison  of  alternative  structures  and  reveals  the  frequency  of  oscilla- 
tion  and  indicates  areas  of  amplitude  sensitivity.  Unless  there  is  a 
known  relationship  between  amplitude  and  the  non-linear  gain  or  amplifi- 
cation factors,  the  amplitude  cannot  be  predicted.   Only  in  cases  where 
there  is  sufficient  low  pass  filtering  is  the  describing  function  an  ade- 
quate relationship  between  gain  and  amplitude.   In  most  oscillator  net- 
works, this  filtering  requirement  is  not  met. 

The  analysis  is  easiest  to  carry  out  if  there  is  one  predominately 
varying  parameter,  i.e.  it  is  a  single  parameter  sensitivity  problem.  A 
multiparameter  sensitivity  index  defined  by  Goldstein  and  Kuo  is  an  ef- 
fective comparison  but  requires  assumptions  concerning  independence  and 
similar  variances  for  all  parameters,  to  be  meaningful. 

Each  oscillator  design  problem  has  its  own  constraints.   In  some 
cases,  minimum  sensitivity  to  parameter  variation  is  desired,  in  another, 
for  a  variable  frequency  oscillator,  a  wide  range  of  frequencies  with 
linear  parameter  variation  is  required.   Consequently,  no  sweeping  gener- 
alities can  be  made  as  to  which  oscillator  circuit  is  best  from  either 
point  of  view.   However,  within  the  context  of  each  design  problem,  there 
is  an  optimum  configuration  for  frequency  sensitivity  to  parameter  varia- 
tions and  such  a  configuration  may  be  determined  by  the  methods  developed, 
provided  that  the  system  characteristic  equations  are  not  of  higher  order 
than  the  third.  Most  of  the  well  known  oscillator  circuits  meet  this 
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criterion.  The  characteristic  equations  of  ten  common  oscillator  cir- 
cuits are  developed  and  the  method  is  easily  adapted  to  variations  of 
these  circuits. 

The  method  employed  suggests  applications  to  the  frequency  charac- 
teristics of  active  and  passive  filters  in  that  the  poles  on  the  joo  axis 
(or  the  points  on  the  £  =  0  contour)  might  equally  well  be  the  poles  in- 
troduced by  a  sinusoidal  driving  function. 

The  worked  example  shows  that  the  single  parameter  sensitivities  are 
functions  of  the  system  parameters.   If  the  sensitivity  with  respect  to 
one  parameter  is  reduced,  it  is  highly  likely  that  the  sensitivity  with 
respect  to  another  has  been  increased.  This  possibility  must  be  consid- 
ered in  any  application  of  the  methods  of  this  paper. 
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APPENDIX  A 

A  BASIC  PROBLEM  ARISING  FROM  THE 

STABILITY  ANALYSIS  OF  LIMIT  CYCLES 

IN  THE  PARAMETER  PLANE 

Introduction 
In  the  complex  frequency  plane,  a  fundamental  criterion  of  absolute 
stability  is  that  there  shall  be  no  roots  of  the  system  characteristic 
equation  lying  in  the  right  half  plane.   The  imaginary  axis,  or  £  a  0, 
line  in  the  complex  plane  is  the  boundary  between  the  regions  of  stabil- 
ity and  instability.   This  boundary  can  be  mapped  into  the  parameter 
plane  as  a  £  a  0  contour.   In  parameter  plane  theory,  this  contour  is 
considered  to  be  the  boundary  between  regions  of  stable  solutions  and 
unstable  regions.   However,  an  example  arising  out  of  the  analysis  of 
the  stability  of  limit  cycles  in  the  oscillator  systems,  in  this  thesis, 
indicates  that  this  is  not  the  case  and  that  there  are  other  boundaries, 
other  than  the  f   =  0  boundary  in  the  parameter  plane  to  mark  off  regions 
of  stability  and  instability. 

The  Particular  Problem 
The  reduced  characteristic  equation  for  an  RC  oscillator  using  a 
bridged  "T"  network  as  a  feedback  network  is  given  by 

(a-ag)p2+  (2f3a-2*e)p+  (l-kx-a£)  =  0  (40) 

The  first  step  is  to  find  the  parametric  equations  of  the  £  =  0  contour 
in  the  a,  (3  parameter  plane.   By  the  method  employed  in  the  thesis,  the 
argument  is  offered  that  £  =  0  for  all  cases  when  the  coefficient  of  the 
first  order  term  in  (40)  vanishes.   Hence  the  non-parametric  equation  of 
the  contour  is  given  by: 
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2  +  2a   -  2c*B  =  0  (68) 

from  which  it  is  seen  that 

P-«  +  f  (69> 

This  is  the  equation  of  a  rectangular  hyperbola  with  asymptotes  a   =  0 
and  B  =  1.5.   Substituting  s  =  j^f  in  the  characteristic  equation  and 
setting  the  real  and  imaginary  parts  equal  to  zero  has  the  same  conse- 
quence. This  is  by  way  of  insisting  that  the  solutions  be  a  pair  of 
purely  imaginary  complex  roots.   This  leads  to  the  equations: 

(a-aB)X2=  (1+a-aB)  (70) 

(2+3a-2aB)  a  0  (68) 

The  parametric  equations  of  the  £  =  0  contour  are  then 

a  =  -2£z  (41) 

ft    l+Z*?  (,o\ 

3  =  ^»?r  (42) 

Lest  there  be  some  argument  that  this  does  not  come  from  the  general 
parameter  plane  method,  the  equations  will  be  developed  from  the  gener- 
alized parameter  plane  method.   Table  2  shows  the  development  of  the 
necessary  intermediate  functions  to  find  the  parametric  equations  of  the 
£  ■  0  contour.   It  is  seen  that  the  parameter  plane  method  in  general 
maps  the  imaginary  axis  into  the  rectangular  hyperbola  developed  before 
from  basic  concepts  (and  a  simpler  method)  and  into  a  point  at  infinity, 

(0,  CO). 

If  the  £  a  0  contour  is  the  boundary  between  regions  of  stability 
and  instability,  a  set  of  test  points  should  reveal  the  existence  of 
stable  and  unstable  solutions  in  accordance  with  this.  This  is  done  for 
5  points,  a  =  10,  B  =  0,  0.5,  1,  1.5,  2.  The  Routh-Hurwitz  stability 
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criterion  is  applied  for  each  set  of  coefficients  determined.  This  de- 
velopment is  shown  in  Table  3. 

Conclusion 
The  £  =  0  curve  does  not  provide  the  only  boundary  in  the  parameter 
plane  between  regions  of  stable  and  unstable  solutions  to  the  system 
differential  equation. 
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Table  2 
Parameter  Plane  Analysis  -  Example 

ak  =  V  +  ^  +  V*  +  dk 
xi  =  ('1)kVkVi(C)   x2  -  (-Dkxk^k\(C) 

k  k 


k                        b 

£ 

h 

d 

2                        1 

0 

-1 

0 

1                        3 

0 

-2 

2 

0                        1 

0 

-1 

1 

2       ^ 
B,   m   U)      -    1 

B2  s   -  3u) 

c1=o 

c2=o 

Dx  s   -   1 

D2  a    -    2lD 

HL  =    -  a)2  +    1 

a 
b 
c 

3 

3    U)       -    U) 

=    0 
=   0 

H2  =   2u) 
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,   -  2w3  - 

U) 

e 

=    f  *   2u)3 

J 

3 

»)    7^     0 

e 

or  -  — 
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-e  ±  ft/e2   ; 
1,2"                 2a 

-  4ac  =   0 

Qt\    w    0 

2u)2 
a2ss    1   -, 

»2 

M  1 , 2  "        H,a . 

2  +   Di 
2+   Cl     " 

B2pfl,2 
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+   D2 
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Pj_  =  °° 

P2=      2u)2 
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Table  3 
Routh  Stability  Criterion  -  Example 

Substitution  of  points  in  parameter  plane  and  application  of 
Routh-Hurwitz  stability  criterion.   For  second  order,  the  Routh  array 
in  general  is 

a2  a0 

*1   ° 

a0  ° 
0 

i.e.  if  there  is  no  sign  change  in  the  sequence  of  coefficients  the 

system  is  stable.   Number  of  sign  changes  is  number  of  roots  in  right 

half  plane,  a  a  10  throughout. 


g. 

fl 

fl 

fo 

Roots  in  RHP 

0 

10 

32 

11 

0  stable 

0.5 

5 

22 

6 

0  stable 

1.0 

0 

12 

1 

transition 

1.0  - 

0  + 

12  + 

1  + 

0  stable 

1.0  + 

0  - 

12  - 

1  - 

1  unstable 

1.5 

-  5 

2 

-  4 

2  unstable 

2.0 

-10 

-  8 

-  9 

0  stable 

Figure  21  shows  the  £  a  0  contour  and  the  points  considered. 
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Figure  21 
£  =  0  Contour  for  Bridged  "T"  Network  Oscillator 
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Suggestions 
The  mapping  of  the  complex  frequency  plane  into  the  parameter  plane 
is  such  that  some  points  map  into  lines.   It  is  possible  for  roots  to 
appear  in  the  right  half  plane  after  leaving  the  left  half  plane  at  in- 
finity, hence,  the  point  at  infinity  must  be  carefully  mapped  in  the 
parameter  plane.  At  the  origin  of  the  complex  frequency  plane  it  is  im- 
possible to  speak  of  the  value  of  £  since  £  is  indeterminate  at  the  ori- 
gin. This  raises  the  possibility  that  the  origin  must  be  separately 
mapped  into  the  parameter  plane,  to  dilineate  areas  of  different  numbers 
of  roots  in  the  left  half  plane. 

Mapping  the  point  at  infinity  yields 
(a-ap)  ■  0 
which  is  the  pair  of  lines  a   =  0  and  f3  a  1.   Both  of  these  lines  are 
significant  boundaries  when  considering  the  Routh  criterion. 

Mapping  the  origin  yields 
(1  +  a   -  ap)  =  0 
This  maps  into  another  rectangular  hyperbola  asymptotic  to  or  a  0  and 
P  ss  1.   This  line  also  yields  a  boundary  between  regions  of  varying  num- 
bers of  roots  in  the  left  half  plane. 

If  the  Jacobian  is  considered  in  accordance  with  p2]  ,  it  is  seen 
that  there  is  a  singularity  at  a   =  0.  This  is  confirmation  of  the  above 
boundaries. 

It  is  evident  that  one  must  be  very  careful  not  to  generalize  about 
the  absolute  nature  of  the  £  is  0  contour  as  the  boundary  between  regions 
of  stability  and  instability.  This  represents  a  boundary  only  for  the 
case  of  a  pair  of  complex  roots  moving  into  the  right  half  plane  of  the 
complex  frequency  plane. 
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